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ABSTRACT 

We present early results from a Hubble Space Telescope (HST) WFC3/IR imaging survey of star- 
forming galaxies in the redshift range 1.5 < z < 3.0. When complete, this survey will consist of 42 
orbits of F160W imaging distributed amongst 10 survey fields on the line of sight to bright background 
QSOs, covering 65 arcmin^ to a depth of 27.9 AB with a PSF FWHM of 0.18". In this contribution, 
we use a subset of these fields to explore the evolution of the galactic stellar mass-radius relation for 
a magnitude-limited sample of 102 spectroscopically-confirmed star forming galaxies ((SFR) ~ 3OM0 
yr^^) with stellar mass ^ IO^^Mq. Although the light profile of these galaxies often has an 
irregular, multi-component morphology, it is typically possible to describe the brightest component 
with a Sersic profile of index n 1. The circularized half-light radius Tq of the brightest component 
is on average {r^) — 1.66 ± 0.79 kpc (i.e., ~ 50 — 70% the size of local late-type galaxies with similar 
stellar mass), consistent with recent theoretical models that incorporate strong feedback from star 
forming regions. The mean half-light radius increases with stellar mass and, at fixed stellar mass. 



evolves with cosmic time as ~ (1 



-1.42 



suggesting that high redshift star forming galaxies may 



evolve onto the local stellar mass-radius relation by redshift z ^ 1. 

Subject headings: galaxies: high-redshift — galaxies: fundamental parameters 



galaxies: structure 



1. INTRODUCTION 

The redshift range z ~ 1.5 — 3.0 is a critical transition 
period in the evolution of galaxies. It is at these redshifts 
when modern-day galaxies are thought to form the ma- 
jority of their stars (e.g., Reddy et al. 2008), fueled by 
typical star formation rates (SFR) ^ 3OM0 yr~^ (Erb et 
al. 2006a) comparable to those seen in local starbursts 
and driving strong enriched winds into the surrounding 
intergalactic medium (see, e.g., Steidel et al. 2010). Pre- 
vious work has shown that this star formation generally 
occurs in galaxies with irregular, often multi-component 
morphologies (see, e.g., Conselice et al. 2005; Ravin- 
dranath et al. 2006; Law et al. 2007; and references 
therein) at 2; > 2, and these galaxies must experience 
a strong morphological transformation in order to form 
the Hubble "tuning fork" ensemble by z ~ 1.0 (e.g., Pa- 
povich et al. 2005; Ravindranath et al. 2004). 

One of the basic properties that is useful for constrain- 
ing models of galaxy formation and stellar feedback pro- 
cesses is the stellar half-light radius To, and the evolu- 
tion of this quantity with stellar mass and cosmic time. 
Driven by the visible-wavelength surveying efficiency of 
the ACS camera on board the Hubble Space Telescope 
{HST), early efforts to characterize the morphologies of 
galaxies at z ~ 1.5 — 3 (e.g., Abraham et al. 1996; Con- 
sehce et al. 2004; Lotz et al. 2006; Papovich et al. 2005; 
Law et al. 2007) have generally focused on rest-frame 
UV emission tracing the regions of active star formation. 
While rest-frame UV and rest-frame optical morphology 
are generally similar for many star forming galaxies (e.g.. 
Bond et al. 2011), there can be significant morphological 
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differences between the evolved stellar population (traced 
by rest-optical emission) and the sites of active star for- 
mation (traced by rest-UV emission) for some high-mass 
galaxies (e.g.. Toft et al. 2005; Law et al. in prep.) 
that complicates the establishment of a robust relation 
between stellar half light radius and stellar mass. 

Since rest-frame optical continuum emission is red- 
shifted into the near-IR for galaxies at z > 1 how- 
ever, it has been challenging to constrain Tc for a sta- 
tistically significant sample of galaxies. Typical z ~ 2 
galaxies are either unresolved or poorly-resolved in see- 
ing limited ground-based imaging, and studies using the 
HST/mCMOS camera (e.g., Conselice et al. 2011) or 
ground-based adaptive optics fed imagers (e.g., Carrasco 
et al. 2010) have too narrow a field of view to permit 
efficient surveys of large numbers of galaxies. Nonethe- 
less, such efforts to characterize the evolution of us- 
ing ground-based instruments and/or the HST/NIGMOS 
camera have been made by (e.g.) Papovich et al. (2005), 
Franx et al. (2008), Toft et al. (2009), van Dokkum et 
al. (2010), and Mosleh et al. (2011), generally finding 
that galaxies at z ~ 2 were significantly smaller at fixed 
stellar mass than in the local universe. 

With the advent of the new WFC3 camera onboard 
HST, it has recently become practical to perform wide- 
field morphological surveys in the near-IR that trace 
rest-frame optical emission from galaxies at z > 1 (e.g., 
Cameron et al. 2010; Cassata et al. 2010). Due to a 
combination of observational limitations and K-hand se- 
lection techniques however, these and previous rest-frame 
optical studies have generally focused upon galaxies with 
stellar mass Af* > 5 x IO^^Mq that are not representa- 
tive of the bulk of the z ~ 2 — 3 star forming galaxy 
population. In addition, many studies have relied prin- 
cipally upon photometric redshifts, which typically have 
large uncertainties (Az/ {l+z) > 0.06; van Dokkum et al. 
2009) at 2 > 1.5, and in individual cases can sometimes 
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fail catastrophically. 

We have therefore undertaken an HST/WFC3 imaging 
survey to map the rest-optical morphology of a large sam- 
ple of optical color selected, spectra scopically confirmed 
galaxies in the redshift range z ~ 1.5 — 3.0 with stellar 
masses in the range ~ 10^ — IQ^^Mq. A complete de- 
scription of this survey and a full analysis of the morpho- 
logical properties of these galaxies will be given by Law 
et al. {in prep.): In this contribution we present prelimi- 
nary results on the mass-radius relation for z ~ 1.5 — 3.0 
star forming galaxies based on early observations taken 
as part of our larger imaging survey. In §2 we present the 
observational data and describe our galaxy sample, out- 
lining our results in §3 and discussing their implications 
for galaxy formation models in Section 4. We assume a 
standard ACDM cosmology in which Hq = 71 km s~^ 
Mpc-\ = 0.27, and I^a = 0.73. 

2. DATA, SAMPLE, AND ANALYSIS 

The HST WFC3/IR camera was used to obtain data 
as part of the Cycle 17 program GO-11694. When com- 
plete, this program will consist of 42 orbits of integra- 
tion in the F160W filter divided amongst 14 pointings 
in 10 different survey fields centered on lines of sight to 
bright background QSOs (z ~ 2.7). Each of our point- 
ings consists of three orbits containing three individual 
exposures of 900 seconds each, for a total of 8100 sec- 
onds per pointing. In this contribution, we present early 
results based on the first 6 pointings (obtained between 
Oct 2009 and Aug 2010) in the Q1009+29, Q1217-f49, 
Q1549-I-19, Q1700+64, and Q2343-I-12 fields, which cover 
a total area of 28 arcmin^. Since these fields are dis- 
tributed widely across the sky we expect cosmic variance 
in our combined sample to be greatly reduced relative to 
comparable surveys over contiguous regions of sky. 

These data were reduced using MultiDrizzle (Koeke- 
moer et al. 2002) with a pixel scale of 0.08 arcsec pixel"^ 
and a pixel droplet fraction of 0.7, resulting in the clean- 
est and narrowest PSF while ensuring that the RMS vari- 
ation of the final weight map was less than ^ 7% across 
the field of view. The weight map produced by Mul- 
tiDrizzle was used to construct an RMS map, scaling 
by a correction factor Fa — 0.3933 (see discussion by 
Casertano et al. 2000) to account for correlation of the 
interpixel noise. We find that the imaging data reach a 
limiting depth of 27.9 AB for a 5cr detection within a 
0.2 arcsec radius aperture, and have a PSF FWHM of 
0.18 ± 0.01 arcsec estimated from unsaturated stars. 

Our galaxy sample is drawn from rest-UV color- 
selected catalogs described by Steidel et al. (2003, 2004) 
and Adelberger et al. (2004). These catalogs were 
constructed based on deep ground-based UnGTZ imag- 
ing data. With sizes of < 1 arcsec the target galaxies 
are effectively unresolved in these ground-based imag- 
ing data, and we therefore do not expect our sample 
to be affected by any significant morphological or sur- 
face brightness biases. We consider in this work only 
those 146 galaxies that have been spectroscopically con- 
firmed (using Keck/LRIS rest-UV spectroscopy) to lie in 
the redshift range z = 1.5 - 3.0, i.e., the "BM", "BX", 
and "LBG" (strictly "C", "D", "M", and "MD") color- 
selection criteria described by Steidel et al. (2003, 2004). 
With an effective wavelength of A = 15369 A, the F160W 
filter approximately traces from B- to y-band rest-frame 



wavelengths over the redshift range of our sample. All 
galaxies that have spectroscopic redshifts are detected 
in the WFC3 imaging data, down to a faint-magnitude 
limit of F160W^ ^ 25.2 AB. Since quantitative morpho- 
logical measurements become unreliable in these data 
fainter than F160W^ ^ 24.0 AB however (see discussion 
by Law et al. in prep.), we impose an apparent magni- 
tude requirement that F160W < 24.0. Additionally, we 
exclude from our sample 6 galaxies which are known to 
contain AGN or faint QSO on the basis of either photo- 
metric or spectroscopic indicators. This combination of 
selection criteria results in a final sample of 102 galaxies 
(see Figure [TJ. 

F160W^ magnitudes for the target galaxies were ob- 
tained using SExtractor (Bertin & Arnouts et al. 1996) 
and corrected for nebular line emission as detailed by 
Law et al. (in prep.)^ In addition to the deep ground- 
based UnGTZ imaging data that forms the backbone of 
the color-selected galaxy catalog, most fields also have 
ground-based J- and/or Ks-hand imaging, and in some 
cases Spitzer IRAC and/or MIPS photometry. These 
data were combined to obtain stellar mass estimates by 
performing a spectral energy distribution (SED) fitting 
analysis using Chariot & Bruzual (2011, in prep.) stel- 
lar population synthesis models with a Chabrier (2003) 
initial mass function (IMF). Stellar masses for the star- 
forming galaxy sample vary from M^, = 10^ — 10^^ Mq, 
although our sample is most well-sampled in the mass 
range = lO^-^ - 10^°-^Mq 

Effective (i.e., half-light) radii were obtained using the 
GALFIT 3.0 program (Peng et al. 2010) to fit a 2D Sersic 
powerlaw profile 



S(r) = Ee exp 




(1) 



to the F160W galaxy morphology, using isolated bright 
(i^l60M^ < 20 AB) but unsaturated stars in each field 
for the PSF model. GALFIT calculates the semi-major 
axis radius a by default, we convert this to a circularized 
effective radius using the formula ro — ay^Jb/a) (where b 
is the semi-minor axis). The PSF FWHM is 0.18 arcsec, 
or about 1.5 kpc at z 2, but since GALFIT decon- 
volves the PSF it is possible to measure effective radii 
smaller than this. Using the method described by Toft 
et al. (2007), we find a 3a resolution limit of 0.073 arcsec, 
corresponding to r-e^so- = 0.62 kpc at z = 2. Nine galax- 
ies are computed to have effective radii smaller than this 
value and are thus consistent with a ^ 9% unresolved 
point source fraction. 

Since many galaxies have an irregular, multi- 
component morphology it is challenging to know how 
to define the characteristic radius of these systems, and 
whether it should describe the size of the brightest com- 
ponent clump or the overall size of the system. In the 
present contribution, in galaxies for which there was 
more than one component to the light profile (i.e., a 
single-component Sersic model resulted in a significant 
residual), we used multiple Sersic components to describe 
the galaxy and adopted the circularized effective radius 

* In brief, we estimate approximate nebular line fluxes from the 
UV-derived star formation rate, and adopt the full amount of the 
correction as our uncertainty. 
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TABLE 1 

Mean Circularized Effective Radius 



re 1.66 ± 0.10 kpc 1.51 ± 0.13 kpc 1.84 ± 0.14 kpc 1.5-3.0 

re 1.90 ±0.18 kpc 1.79 ± 0.20 kpc 2.10 ±0.34 kpc 1.5-2.0 

re 1.65 ±0.13 kpc 1.45 ± 0.20 kpc 1.80 ±0.17 kpc 2.0-2.5 

re 1.22 ±0.12 kpc 1.01 ±0.14 kpc 1.50 ±0.13 kpc 2.5-3.0 

f-e/rsDSS 0.61 ±0.04 0.62 ±0.05 0.61 ± 0.05 1.5-3.0 

ro/rsDSS 0.73 ±0.07 0.74 ±0.09 0.72 ±0.12 1.5-2.0 

rc/rsDSS 0.59 ±0.04 0.59 ± 0.08 0.58 ± 0.05 2.0-2.5 

r-o/rsDSS 0.45 ±0.04 0.41 ± 0.05 0.51 ± 0.04 2.5-3.0 



of the primary component (i.e., the brightest in F160W 
flux) as the of the galaxy. The primary component is 
well-defined for most of our galaxies. In ^ 10% of cases, 
however, the primary and secondary components have 
comparable magnitudes, and our stellar masses (which 
were determined from photometry with angular resolu- 
tion too poor to distinguish individual clumps) are likely 
overestimates of the mass of the brightest clump. As dis- 
cussed in fJH however, our results are largely insensitive 
to this complication. 

3. RESULTS 

As illustrated in Figure [U galaxies range in effective 
radius from unresolved to = 5.5 kpc, with a Sersic 
index n < 2.5 in all but three cases (two of which had Tc 
values consistent with an unresolved point source). Typi- 
cal stochastic uncertainties in and n vary as a function 
of magnitude and morphological type, but are generally 
~ 3% and 8% respectively. We caution, however, that 
while the light profile of the primary component of these 
galaxies is therefore "disklike" in the sense that the light 
profile falls off more slowly with radius than for bulge- like 
profiles {n — A), these are generally not classical disks in 
a morphological or kinematic (e.g., Forster-Schreiber et 
al. 2009; Law et al. 2009) sense. 

Plotting the circularized effective radius Vc versus stel- 
lar mass in Figure [5J we observe that although there is a 
large spread in To for individual galaxies, a mass-radius 
relation is generally in place for star-forming galaxies at 
^ lO^^M© since at least z ~ 3.0. Binning our sample 
by redshift and stellar mass (see Table [Ij , we calculate 
that (r-o) = 1.01±0.14 (1.50±0.13) kpc for galaxies in the 
mass range = 10^-^-^° °{10^°-°-^^-^)Mq respectively 
at redshift z = 2.5 — 3.0, increasing with cosmic time to 
(re) = 1.45 ± 0.20 (1.80 ± 0.17) kpc by z = 2.0 - 2.5, and 
to (re) = 1.79 ± 0.20 (2.10 ± 0.34) kpc by z = 1.5 - 2.0. 

4. DISCUSSION 

The radii that we find for galaxies in the redshift 
z ~ 2 — 3 universe (re ~ 1 — 2 kpc for galaxies with stel- 
lar masses M, ~ 10^°Mq) are generally consistent with 
the theoretical predictions of Sales et al. (2010), and 
favor their "WF2Dec" model in which relatively strong 
feedback from star forming regions results in the efficient 
removal of gas from galaxies via an outfiowing wind with 
velocity ~ 600 km s~^ . Such peak outflow velocities are 
generally consistent with observations (see, e.g., Steidel 
et al. 2010). Dividing the radii of each of our galaxies 
by the local value at corresponding stellar mass based 
on the late-type (i.e., n < 2.5) relation found by Shen et 
al. (2003) for ~ 140000 galaxies in the Sloan Digital Sky 



Survey (SDSS), we estimate (see Table [T]) that the radii 
of z ~ 2—3 galaxies are on average (^e/rsDSs) ^ 50—70% 
the size of comparable-mass galaxies in the low redshift 
universe. This is somewhat larger than previous esti- 
mates for high-mass (M* > 5 x lO^^M©) star forming 
and quiescent galaxy populations. Toft et al. (2009) 
for instance found that for a sample of 225 galaxies 
with photometric redshifts from the FIREWORKS cat- 
alog, both star-forming and quiescent galaxies at z ~ 2 
were significantly smaller at fixed stellar mass than in 
the local universe, with (re/rsDSs) = 0.51 ± 0.02 and 
(?'c/''SDSs) = 0.34 ± 0.02 respectively. Similarly, Franx 
et al. (2008) found {rJrsDSs) = 0.52 ± 0.06 for a K- 
selected sample of z ~ 2 galaxies while van Dokkum et 
al. (2008) found (re/rsoss) = 0.17 for a K-selected sam- 
ple of galaxies without nebular emission lines. 

Toft et al. (2009) postulated that the proximity of 
z ~ 2 galaxies to the local mass-radius relation may be 
a function of the specific star formation rate (sSFR) of 
a given galaxy, with their Figure 3 demonstrating the 
increase in (ro/rsDSs) from ~ 0.1 to ~ 1.0 as sSFR 
increases from < 0.2 Gyr^^ (for quiescent galaxies) to 
~ 1.5 Gyr^^ (for star-forming galaxies). Although con- 
cerns regarding the evolving mass-to-light ratio as a func- 
tion of stellar age tend to complicate such trends, we con- 
firm this statement in the sense that our overall popula- 
tion has both a large sSFR (~ 2 Gyr~^) and a relatively 
large mean (ro/rsDSs) — 0.61 ± 0.04, but note that no 
such relation is present within the range of sSFR probed 
by our star forming galaxy sample. 

The effective radius (as a function of stellar mass) of 
the galaxy population evolves significantly closer to the 
local relation rsoss with decreasing redshift (see Figure 
[3]). Clustering statistics suggest that we are seeing the 
same population of galaxies at all three redshifts (Adel- 
berger et al. 2005), and the size of these galaxies in- 
creases from (r-e/rsDSs) — 0.45 ± 0.04 in the redshift 
range z = 2.5 - 3.0 to 0.59 ± 0.04 by z = 2.0 - 2.5, and 
0.73 ± 0.07 by z = 1.5 — 2.0. Assuming growth of the 
form Te ^ (1 + z)~" between z — 3.0 and z — 1.5, a 
least-squares analysis gives best-fit a = 1.42 ±0.50 (solid 
line in Figure [3]), consistent with the re ^ (1 + z)^^-^ and 
re ^ (1 + z)~^^^^ evolution found by van Dokkum et al. 
(2010) and Mosleh et al. (2011) respectively for massive 
galaxies. While such an evolution is broadly consistent 
with classical theories for disk galaxy evolution (e.g.. Mo 
et al. 1998; see also discussion by Papovich et al. 2005, 
Franx et al. 2008), it suggests that the z ~ 2 — 3 star- 
forming galaxy population might be expected to evolve 
onto the local late- type galaxy relation by z ^ 1, in 



4 



The Mass- Radius Relation at 2: ~ 1.5 — 3.0 




Effective Radius (kpc) Sersie Index n 

Fig. 1. — Histogram of masses log(M/M0), redshifts 2, effective radii rc, and Sersic index n for the 102 galaxies in the z = 1.5 — 3.0 star 
forming galaxy population. The blue/black/red histograms respectively indicate galaxies in the redshift ranges z = 1.5 — 2.0, z = 2.0 — 2.5, 
and z = 2.5 — 3.0. The green histogram represents galaxies of all redshifts in the sample. 
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Fig. 2. — Effective circularized radius function of stellar 

mass M, for the z = 1.5 - 2.0 (blue points), z = 2.0 - 2.5 (black 
points), and z = 2.5 — 3.0 (red points) star forming galaxy samples. 
Crosses, open boxes, and open circles represent galaxies selected 
according to the "BX" , "BM" , and "LBG" color-selection criteria 
respectively, upper limits for unresolved sources are denoted with 
arrows. The filled circles and error bars represent the mean value 
and associated uncertainty in the mean for galaxies in each red- 
shift sample with stellar masses in the range M, = 10^'*~^'''''Mq 
and 1O^''"~1''-^M0. The solid black line indicates the low-redshift 
relation for late-type galaxies from Shen et al. (2003). 

agreement with Barden et al. (2005), who found little 
evidence for evolution of the stellar mass-radius relation 
from z = 1 to the present day. We note that our re- 
sults are robust to our choice to adopt Vc of the brightest 
component as the effective radius for the galaxy while 
stellar masses are derived from the total integrated light 
from all clumps within a given galaxy. If we instead scale 
the stellar mass of the galaxy by the fraction of F160W 
fl-ux contained within the primary component (or indeed, 
if we ignore the comparable-magnitude multi-clump sys- 
tems entirely) our derived value of (^o/rgDss) changes by 
less than 1%, and the index a of the redshift evolution 
changes by < 0.5cr. 

Such a trend might, however, be telling us less about 
the growth of galaxies within their dark matter halos 
than about the evolution of the sites of star formation 
in the young universe. It is not necessarily surprising 
that z ^ 2 star forming galaxies do not adhere precisely 
to the local mass-radius relation for late-type galaxies 
since (despite a similar radial index of their rest-frame 
optical light distribution) they represent drastically dif- 
ferent physical systems. At fixed stellar mass, z ^ 2 — 3 
galaxies tend to have much higher gas fractions (~ 50%, 
Erb et al. 2006b; Daddi et al. 2010) than late-type 
disk galaxies in the nearby universe (< 20%; Leroy et al. 
2008), fueling both high star formation rates ^ 30Mq 
yr~^ and star formation surface densities, and driving 
strong gaseous outflows characteristic of local starburst 
galaxies. Similarly, galaxies at z 2 — 3 are significantly 
more morphologically irregular than typical local galax- 
ies (e.g., Abraham et al. 1996; Conselice et al. 2005; 
Law et al. 2007) and have much higher gas-phase veloc- 
ity dispersions a ^ 70 km (e.g., Forster-Schreiber et 
al. 2009; Law et al. 2009). As z decreases from z ~ 2 to 
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Fig. 3. — Effective circularized radius fraction of the local 

relation (rgDSs) at a given stellar mass as a function of redshift. 
Symbols are as in Figure [2] filled circles represent the mean and 
associated uncertainty of all galaxies in each of the three redshift 
ranges. The solid black line indicates growth of the form To ~ 
(1 -I- z)~^'^^ , as derived from a least-squares fit to the data. 

z ~ 1, gas fractions decrease (e.g., Tacconi et al. 2010), 
morphological properties (e.g., Papovich et al. 2005) and 
kinematics (Wright et al. 2009) begin to resemble the 
z ^ universe, and it is perhaps unsurprising that we 
observe evidence for the typical effective radii of galaxies 
with active star formation beginning to resemble their 
local values as well. 

Indeed, galaxies with properties similar to the z 2 
star forming galaxy sample are not unknown in the low- 
redshift universe, but are simply much more rare. The 
sample of local (z < 0.3) gas-rich supercompact UV lu- 
minous galaxies (ScUVLGs, also sometimes known as 
LBAs) described by Heckman et al. (2005) for instance 
are known to have star formation rates (Hoopes et al. 
2007), morphologies (Overzier et al. 2008), and kinemat- 
ics (Gongalves et al. 2010) similar to z ^ 2 star forming 
galaxies, and it is perhaps unsurprising that they ap- 
pear to have similar rest-optical effective radii as well 
(Overzier et al. 2010). 

Although the morphological sample of galaxies pre- 
sented here is larger than any previous spectroscopically 
confirmed z ^ 2 — 3 star forming galaxy sample with 
mass ^ — IO^^Mq, we caution in closing that it is 
still relatively small and therefore challenging to explore 
robustly the contribution of all possible systematic effects 
to the measured Tq. Such systematics include, but are not 
limited to, bandshifting throughout the redshift inter- 
val, biases due to multi-component morphologies, biases 
arising from differences in magnitude between galaxies 
at different redshifts, etc. These and other topics will 
be addressed in detail by Law et al. {in prep.) using 
the larger sample of 400 galaxies in our full imaging 
survey. 
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